Search for Majorana fermions in multiband semiconducting nanowires 



Roman M. Lutchyn 1 ' 3 , Tudor D. Stanescu 1 ' 2 , and S. Das Sarma 1 
1 Joint Quantum Institute and Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111, USA 
2 Department of Physics, West Virginia University, Morgantown, WV 26506, USA 
3 Microsoft Research, Station Q, Elings Hall, University of California, Santa Barbara, CA 93106, USA 

(Dated: compiled February 22, 2011) 

We study multiband semiconducting nanowires proximity-coupled with an s-wave superconductor. 
We show that when odd number of subbands are occupied the system realizes non-trivial topological 
state supporting Majorana modes localized at the ends. We study the topological quantum phase 
transition in this system and analytically calculate the phase diagram as a function of the chemical 
potential and magnetic field. Our key rinding is that multiband occupancy not only lifts the stringent 
constraint of one-dimensionality but also allows to have higher carrier density in the nanowire and as 
such multisubband nanowires are better-suited for observing the Majorana particle. We study the 
robustness of the topological phase by including the effects of the short- and long-range disorder. We 
show that in the limit of strong interband mixing there is an optimal regime in the phase diagram 
("sweet spot") where the topological state is to a large extent insensitive to the presence of disorder. 

PACS numbers: 03.67.Lx, 71.10.Pm, 74.45. +c 



Looking for the elusive Majorana particles is one of 
the most active and exciting current topics in all of 
physics pp. Although originally proposed as a model for 
neutrinos, the current search for Majorana particles is 
mostly taking place in condensed matter or atomic sys- 
tems [2j [3] where these mysterious particles, which are 
their own anti-particles, emerge as effective quasiparti- 
cles from an underlying fermionic Hamiltonian. Quite 
apart from the intrinsic interest associated with the ex- 
otic Majorana particles, the possibility that they can be 
used in carrying out fault-tolerant topological quantum 
computation [4 by suitably exploiting their non-Abelian 
braiding statistics gives an additional technological im- 
petus in the subject. It has been known for a while [5H5] 
that, under suitable conditions, Majorana particles could 
exist at the ends of ID nanowires in the presence of the 
appropriate superconducting (SC) pairing. Also, it has 
been recently shown that the network of Majorana wires 
can be used for braiding [9] and topological quantum 
computation [10 j. Although the semiconducting (SM) 
nanowires [7| [8] are promising candidates for observing 
the Majorana, experimental realization of these proposals 
is challenging because obtaining strictly ID nanowires is 
a very demanding materials problem [IT] . In this Letter 
we establish that one dimensionality, i.e. the occupancy 
of one only subband in the nanowire, is completely un- 
necessary, and Majorana particles can exist under rather 
general and robust conditions even when several sub- 
bands are occupied in the nanowire. More importantly, 
we prove the remarkable counter-intuitive result that the 
multisubband system is, in fact, better-suited in observ- 
ing the Majorana than the strict ID limit. We carry 
out an analytic theory establishing our main results and 
provide support for it by independent numerical calcu- 
lations. We also study the robustness of the topological 
phase against short- and long-range disorder and show 



that there is an optimal parameter regime where the sys- 
tem is most stable with respect to disorder. We believe 
that our results would go a long way in providing the 
most suitable solid-state system for the eventual obser- 
vation of the Majorana particles. 

In this Letter we propose to study Majorana physics 
in SM quantum well based on, for example, InAs-AlSb 
heterostructure [13] . The active system consists of a SM 
with strong spin-orbit interaction proximity-coupled with 
an s-wave SC, see Fig. [T^l. The rectangular quantum well 
has the dimensions L z , L y and L x as shown in Fig [T^l. 
We consider here the case of a strong confinement in the 
z direction such that L z <C L yi L x so that only the low- 
est subband with respect to the i-axis eigenstates is oc- 
cupied. Then, the single-particle Hamiltonian takes the 
usual form for the 2D SM in the presence of the spin-orbit 
Rashba interaction (h= 1): 

Hsm = J dxdy^l(x,y)H aa ^ a >(x,y) (1) 
d 2 x + d 2 

H = - * V ~ V>~ ia(°xd y - cF y d x ) + V x a x , (2) 

where m*, a and /i are the effective mass, the strength 
of spin-orbit interaction and chemical potential, respec- 
tively. The latter can be controlled using the gate elec- 
trodes [TTJ [13] . The last term in Eq. ([2| corresponds to 
the Zeeman term due to the applied external magnetic 
field aligned along the x-axis, V x = gsM/^BB x /2. Note 
that magnetic field is essential here - it opens up a gap in 
the spectrum at p x = and allows one to avoid fermion 
doubling which is detrimental for the existence of Majo- 
rana fermions. Because of the large g-factor in the SM 
#inAs i$ 35, fairly small in-plane magnetic field B x < IT 
opens up a sizable gap in the spectrum without signifi- 
cantly disturbing the SC. 

We now include the size quantization along ^-direction 
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FIG. 1. (Color online) (a) Schematic plot of the quasi- ID 
nanowire proximity- coupled with an s-wave SC. The rect- 
angular quantum well has the dimensions L z , L y and L x : 
L z <^ L y <^ L x . The nanowire can be top gated to control 
chemical potential in it. The method for fabricating the pro- 
posed quantum well heterostructure based on InAs has been 
demonstrated, see, e.g., Ref. [I2j[l3].(b) Schematic plot of the 
lowest energy subbands due to the transverse confinement. 



assuming that L y <C L x . One can notice that Hamilto- 
nian Q is separable in x — y coordinates and the field 
operator can be written as 



Px ,n y =l,2,... 
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where a r , 



is electron annihilation operator in a state 



n y having momentum p x . The physical parameter regime 
we consider here corresponds to the confinement energy 
along ^/-direction being larger than all the relevant energy 
scales of the Hamiltonian Q so that there are only few 
lowest subbands occupied, see Fig. [TJd. This assumption 
actually corresponds to the typical experimental situa- 
tion in InAs nanowires [13] . For p < 2£? s b, where E^ = 
37r 2 /2m* Ly is the subband energy difference, one can 
project the wavefunction to the lowest two subbands n y = 
1,2 in ([3| and simplify the Hamiltonian 0. By intro- 
ducing the spin-band spinors = (c Px ^, c Px ^ 1 d Px ^, d Px ^) 
where the annihilation operators c Px and d Px correspond 
to a p x ,ny=i and a Px , ny =2, respectively, the single-body 
Hamiltonian becomes H = ^2 Px & (p x )H re a<$>(p x ) with 
i^ re d being defined as 



P 1 — Pz 

H Te d=-^-^-Oia y p x -\-E sh — E hm a x p y + V x a x . (4) 



Here Pauli matrices &i and pi act on the spin and 
band degrees of freedom. The band mixing energy 
^bm corresponds to the expectation value p y opera- 
tor between different band eigenstates, i.e. ^bm = 
/ L "d^ S in(^Hsin(^) = ^. 

We now study topological properties in this regime 
with low number of subbands occupied. We investigate 
here whether Major ana fermions survive and are robust 
in this quasi- ID geometry. The multiband proximity- 



induced SC can be described as 

H SC = [Ancl t cl pxl + A 22 </_^ (5) 
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where the induced SC pairing potentials A^- de- 
pend on the microscopic details of the interface be- 
tween SM and SC, e.g. rough or smooth inter- 
face. In the former case the magnitude of A12 can 
be a sizable fraction of An. Taking into account 
the total Hamiltonian H tot = H re( ± + Hsc we can 
now define the Nambu spinor as follows: ^(p) = 

(c Px t,c Px i,d p ^,d pU ,cl pxV cl pxl ,dl pA ,dl pxl ) . In 
this convention for Nambu spinors the BdG Hamiltonian 
for two subband model reads 

-EbmO-xpy+icTy [p x I A12 1 + A + +p z AJ\(ir y cos <p+ir x sin (p) , 

where Pauli matrices c^, pi and n act on spin, band 
and Nambu degrees of freedom of the spinor ^(p), re- 
spectively; A± = (| An I i IA22D/2 and ip is the SC 
phase. The particle- hole symmetry for i^BdG (|6f reads 
0^BdG(p)0 -1 = — ^BdG(— p), where 6 is an anti- 
unitary operator Q = r x K with K denoting the complex 
conjugation. 

The presence of Majorana modes in the system and the 
corresponding phase diagram can be obtained using topo- 
logical arguments due to Kitaev [5] . Following Ref. [5] we 
introduce Z 2 topological index M (Majorana number): 



M = sgn[Pf5(0)]sgn[Pf£(7r/a)] = ±1, 



(7) 



where ±1 corresponds to topologically trivial/non-trivial 
states. Here antisymmetric matrix B defines the Hamil- 
tonian of the system in the Majorana basis [5]. Rather 
than computing the transformation matrix to the Majo- 
rana basis as was originally done in Ref. [5], we note fol- 
lowing Refs. [3 [15] that the antisymmetric matrix B(p x ) 
can be simply constructed by the virtue of the particle- 
hole symmetry. Indeed, the matrix B(P) = H^&q{P)t x 
needs to be calculated at the particle-hole invariant 
points where i^BdG(^) =^BdG( — P) and B is antisym- 
metric B T (P) = —B(P). In ID there are two such points: 
P = 0. - with - being the momentum at the end of the 

'a a & 

Brillouin zone and a being the lattice spacing. (For the 
continuum model considered here ir/a^oo.) The func- 
tion Pf in Eq. ^ denotes Pfaffian of the antisymmetric 
matrix B. The computation of Pfaffian at P = 7r/a— ^ 00 
is straightforward yielding sgn[Pf£?(7r/a)] = +1. Thus, 
the phase boundary for the transition between topolog- 
ically trivial and non-trivial phases is given by the sign 
change of Pf£?(0) which can only happen when the bulk 
quasiparticle gap becomes zero, i.e Det i^BdG^) =0, see 
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FIG. 2. (Color online) (a) Phase diagram for the two band 
nanowire model as a function of the chemical potential jl and 
external magnetic field V x . The width of the topological re- 
gion is largest at the "sweet" spot: W^60. Here tilde denotes 
re-scaled energy E — E/m*a 2 and A12 = 4. The light and 
dark regions correspond to topologically trivial/non-trivial 
phases, (b) and (c) Quasiparticle excitation gap obtained 
using Eq. |6]) as a function of V x and jl. The solid (red) and 
dashed (blue) lines correspond to A12 =4 and A12 = 1, respec- 
tively. The closing of the gap for A12 =4 (solid red line) is con- 
sistent with the phase diagram shown in (a). The quasipar- 
ticle excitation gap at the "sweet spot" strongly depends on 
the magnitude of A12. We assume here m* = 0.04ra e with m e 
being electron mass and a — O.leVA yielding m*a 2 « 0.6K. 
We used realistic parameters L y = 130nm and E^ = 30, 
E hm — 5 and An = A22 — 4. 



Fig. [2] This is a generic phenomenon since the topologi- 
cal reconstruction of the fermionic spectrum cannot occur 
adiabatically and requires the nullification of the bulk ex- 
citation gap [TTl [18] . For a two-band model Pf£?(0) can 
be calculated analytically 



Pf J B(0) = (K c 2 -i3L + |Ai 2 | 2 +A 2 _-A 2 + - J B sbM + M 2 
— Vx {El-4E shf i+4 (A 2 2 + A 2 _+ M 2 )) 
+ (£ sb (A_+A + )-2A+^) 2 



(8) 



allowing one to compute M. as a function of the physi- 
cal parameters. The phase diagram showing a sequence 
of topological phase transitions for the two subband 
nanowire is shown in Fig. [2^i. We now analyze the phase 
diagram in various regimes. In the limit /i, |V^| <C E s ^ 
we find that PfB(0) ~ — V^ 2 +A^+/i 2 recovering the pre- 
vious results obtained for the single band [TJ [8] . When 
\V X \ <C/i~£ sb we find that PfB(0) « -V*+A 2 22 +(E sh -Li) 2 . 
Thus, the system supports Majorana modes as long as 
I V x I > Y // A22 + (^sb — m) 2 - These results can be intuitively 
understood within weak-coupling approximation since in 
both cases the Fermi level crosses odd number of bands 
in the interval (0, J). The most interesting parameter 
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FIG. 3. (Color online) (a) and (b) Energy spectrum E m for a 
finite-size nanowire obtained by numerical diagonalization of 
Htot = Hsm+Hsc for p, = 15, V x = 15 and V x «8, respectively. 
Here E = E/m*a 2 and m labels eigenvalues of if tot- Inset: 
lowest lying energy states. Majorana zero-energy modes are 
present in (a) and disappear in (b). (c) and (d) Energy spec- 
tra for a given disorder realization with impurity potentials 
shown at the bottom of the panels (c) and (d) corresponding 
to A = 16nm and A = 260nm, respectively. The dark (light) 
colors denote positive (negative) Uoj . c) The spectrum at the 
"sweet spot" for the short-range disorder shown at the bottom 
(Uo = 100). Inset: lowest-lying eigenstates. The topological 
phase is robust against short-range disorder, i.e. the disorder 
affects extended states but leaves Majorana modes intact, d) 
The spectrum at the "sweet spot" with long-range disorder 
for [To = 100 (squares) and Uo = 25 (diamonds). The impurity 
potential t/i mp is shown at the bottom of the panel. Inset: 
lowest-lying eigenstates for Uo = 100. The topological phase 
is stable as long as Uo < W/2. For Uo = 100 additional Majo- 
rana modes are localized at the impurities and the topological 
phase collapses by splitting into fragments of topological and 
non-topological regions, see inset. Both types of disorder lead 
to the emergence of the additional subgap states localized at 
the ends. Here we used parameters specified in Fig. [2] 



regime is fi ~ which corresponds to the "sweet 

spot" in the phase diagram, see Fig. [2^i. At this point the 
system is to a large extent insensitive to chemical poten- 
tial fluctuations and, thus, this regime provides a promis- 
ing route to realizing a robust topological SC phase. At 
H = E s }j/2 the width of the topologically non-trivial re- 
gion is given by E sh /2 - A 12 < \V X \ < E sh /2 + A 12 to 
a leading order in 1/E^. This is non-perturbative re- 
sult and the SC state emerging here is determined by the 
strong inter band mixing due to A 12. The presence of a 
sizable A12 is crucial for the topological stability of the 
non-trivial SC phase and the magnitude of the quasipar- 
ticle excitation gap at the "sweet spot" strongly depends 
on the value of A 12, see Figs. [2)3 and c. 

In order to establish the robustness of the topologically 
non-trivial phase near the "sweet spot" we have done in- 
dependent numerical simulations for a finite multiband 
nanowire with L x ^10liir and L y ~ 0.1/im. The results 
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obtained by numerical diagonalization of the real-space 
Hamiltonian H tot =#sm+#sc are shown in Fig.JS^i and 
b. One can notice that at the "sweet spot" (V x « 15 
and jl w 15) there is a pair of Major ana zero-energy 
states whereas for a smaller magnetic field (V x ~ 8) cor- 
responding to the trivial phase the zero energy modes 
disappear corroborating the phase diagram shown in 
Fig. [2^l. Furthermore, at the "sweet spot" the zero en- 
ergy states are well-separated from the continuum. In- 
deed, as shown in Fig. [3^i the minigap E mn constitutes a 
sizeable fraction of the induced SC gap, E mn ~lK. Thus, 
Major ana modes in quasi- ID nanowires are very robust 
against thermal fluctuations which makes these systems 
very advantageous for the topological quantum compu- 
tation. We also studied the robustness of the topological 
phase against disorder by adding the impurity potential 
^imp(y) = J2j u 0j e *\^r^\/d X) to the Hamiltonian H tot . 
Here r is a vector in x—y plane, rj are random positions of 
the impurities, A is the screening length, Uqj = ±[/o is the 
impurity potential with random sign but constant mag- 
nitude /To, and d is the cutoff regularizing 1/r potential 
at short distances. We considered here two types of dis- 
order mimicking short-range impurities (A = 16nm) and 
long-range chemical potential fluctuations (A = 260nm), 
see Fig. |3]^,d. In the former case, the topological phase 
is very robust against disorder even if \Uo\ ^> ^ s b> see 
Fig. This can be qualitatively understood as follows: 
at a given position the disorder potential leads to a for- 
mation of two Majorana modes localized at the impurity. 
Because these Majorana states are close to each other 
they hybridize and form conventional subgap states and 
do not affect Majorana modes at the ends even if the 
impurity is fairly close to the edge, see Fig. 3c. On the 
other hand, the long-range disorder is more dangerous. 
In Fig. |3]i we show energy spectrum for two cases: Uq 
smaller (U = 25) and larger (U = 100) than W/2, see 
Fig. [2|i. For Uo < W/2, the topological phase is sta- 
ble, i.e. the disorder can suppress excitation gap but 
does not affect Majorana modes. On the other hand, if 
Uq>W/2, the disorder effectively creates inhomogeneous 
wire with many topological and non- topological regions, 
see Fig. i.e. each topological segment now becomes 
much smaller allowing for strong mixing of the Majorana 
modes at the opposite ends. Thus, our simulations ex- 
plicitly demonstrate the importance of working at the 
"sweet spot" where the width of the topological region 
W is maximized and the topological phase is most robust 
against long-range disorder. 

To conclude, we have derived the topological phase 
diagram for the existence of Majorana particles in a re- 
alistic quasi- ID semiconductor wire in the presence of 
multisubband occupancy. Unexpectedly, we find robust 
and experimentally feasible "sweet spots" in the chemical 
potential- Zeeman splitting phase diagram where Majo- 



rana modes should stabilize at the ends of the wire. The 
great advantages of our proposed structure in detecting 
Majorana particles are (i) its materials flexibility (i.e. no 
need to impose one dimensionality or single channel con- 
straint), and (ii) its immunity to density (or chemical 
potential) fluctuations and disorder. The calculation of 
the energy spectrum for realistic experimental settings 
suggests the possibility to test our theoretical results us- 
ing local tunneling experiments, see Ref. [16]. Tunneling 
of electrons to the ends of the nanowire would reveal a 
pronounced zero-bias peak when the system is in topo- 
logically non-trivial phase. This zero bias peak will dis- 
appear in the trivial phase. 

Note added. While finishing this manuscript we be- 
came aware of Refs. [3j [19] where multichannel general- 
ization of the spinless p-wave SC state was studied. 

This work is supported by DARPA-QuEST and JQI- 
NSF-PFC. 
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